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1 Introduction 

The stochastic calculus with respect to the fractional Brownian motion (fBm) has now a 
long enough history. Since the nineties, many authors used different approaches to develop a 
stochastic integration theory with respect to this process. We refer, among of course many 
others, to [T], [8], |12j or [9]. The reason for this tremendous interest in the stochastic anal- 
ysis of the fBm comes from its large amount of applications in practical phenomena such as 
telecommunications, hydrology or economics. 

Nevertheless, even fBm has its limits in modeling certain phenomena. Therefore, several 
authors introduced recently some generalizations of the fBm which are supposed to fit better 
in concrete situations. For example, we mention the multifractional Brownian motion (see e.g. 



[2]), the subfractional Brownian motion (see e.g. [5]) or the multiscale fractional Brownian 
motion (see [3]). 

Here our main interest consists in the study of the bifractional Brownian motion 
(bifBm). The bifBm has been introduced by Houdre and Villa in [13J and a stochastic analysis 
for it can be found in [19]. Other papers treated different aspects of this stochastic process, 
like sample pats properties, extension of the parameters or statistical applications (see [6], [I], 
[22j or [IT]. Recall that the bifBm B ' is a centered Gaussian process, starting from zero, 
with covariance function 

R H ' K (t, s) := R(t, s) = ± + s™) K - \t - s? HK ) (1) 

where the parameters H, K are such that H 6 (0,1) and K £ (0,1]. In the case K = 1 we 
retrieve the fractional Brownian motion while the case K = 1 and H = | corresponds to the 
standard Brownian motion. 

The process B H)K is ifif-selfsimilar but it has no stationary increments. It has Holder 
continuous paths of order 5 < HK and its paths are not differentiable. An interesting property 
of it is the fact that its quadratic variation in the case 2HK = 1 is similar to that of the standard 
Brownian motion, i.e. [B ' ]t = est. x t and therefore especially this case {2HK = 1) is very 
interesting from the stochastic calculus point of view. 

In this paper, our purpose is to study multidimensional bifractional Brownian motion 
and to prove Ito and Tanaka formulas. We start with the one dimensional bifBm and we 
first derive an Ito and an Tanaka formula for it when 2HK > 1. We mention that the ltd 
formula has been already proved by [15J but here we propose an alternative proof based on 
the Taylor expansion which appears to be also useful in the multidimensional settings. The 
Tanaka formula is obtained from the Ito formula by a limit argument and it involves the 
so-called weighted local time extending the result in [7J. In the multidimensional case we first 
derive an Ito formula for 2HK > 1 and we extend it to Tanaka by following an idea by Uemura 
|20j . |21j ; that is, since \x\ is twice the kernel of the one-dimensional Newtonian potential, i.e. 
^A|x| is equal to the delta Dirac function 5(x), we will chose the function U(z), z G M. d which is 
twice of the kernel of d-dimensional Newtonian (or logarithmic if d = 2) potential to replace \x\ 
in the cf-dimensional case. See the last section for the definition of the function U. Our method 
is based on the Wiener-Ito chaotic expansion into multiple stochastic integrals following ideas 
from [14J or [10] . The multidimensional Tanaka formula also involves a generalized local time. 
We note that the terms appearing in our Tanaka formula when d > 2 are not random variables 
and they are understood as distributions in the Watanabe spaces. 



2 Preliminaries: Deterministic spaces associated and Malliavin 
calculus 

Let ^B^' K ,t £ [0,T]^ be a bifractional Brownian motion on the probability space {Q.,J-,P). 

Being a Gaussian process, it is possible to construct a stochastic calculus of variations 
with respect to B ' . We refer to [l], [17] for a complete description of stochastic calculus 
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with respect to Gaussian processes. Here we recall only the basic elements of this theory. 

The basic ingredient is the canonical Hilbert space Ti associated to the bifractional 
Brownian motion (bifBm). This space is defined as the completion of the linear space £ 
generated by the indicator functions l[o,t]>£ S [0,T] with respect to the inner product 

(l[0,t]> \o,s])h = R(t, s). 

The application ip G £ — > B{ip) is an isometry from £ to the Gaussian space generated by 
B H ' K and it can be extended to H. 

Let us denote by S the set of smooth functionals of the form 

F = f(B(<p 1 ),...,B(<p n )) 

where / G C£°(M n ) and (fi G 7i. The Malliavin derivative of a functional F as above is given 
by 

D B *'*F = J£(i%i), • • • , BWM 

i=l 1 

and this operator can be extended to the closure B m ' 2 {m > 1) of 5 with respect to the norm 

\\F\\l, 2 := E\F\* + E\\D BH - K Ff H + ... + E\\D bH ' K > ™F||^ m 

where 7^® m denotes the m fold symmetric tensor product and the mth derivative D bH K ,m is 
defined by iteration. 

The divergence integral 5 B ' is the adjoint operator of D B ' ^ . Concretely, a random 
variable u G L 2 {Q;TL) belongs to the domain of the divergence operator (Dom(5 B ' )) if 



E 



(D bH K F, u) 



H 



< cII^IIl 2 ^) 



for every F G S. In this case 5 bH K (u) is given by the duality relationship 

E(F6 BH ' K (u)) = E(D BH ' K F,u) n 

for any F G B 1 - 2 . It holds that 

E5 BH ' K {uf = E\\uf n + E(D BH ' K u, (D BH ' K uy) n ® H (2) 

where (D bH ' K u)* is the adjoint of D bHK u in the Hilbert space Ti^Ti. 

Sometimes working with the space 7i is not convenient; once, because this space may 
contain also distributions (as, e.g. in the case K = 1, see [18]) and twice, because the norm in 
this space is not always tractable. We will use the subspace \H\ of Ti which is defined as the 
set of measurable function / on [0, T] with 

f H] := £ £ \f(u)\ \f(v)\ ^-( U ,v)dudv < oo. (3) 
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It follows actually from [15] that the space \7i\ is a Banach space for the norm || • and it 
is included in TC. In fact, 



L 2 ([0,T]) C Lttk([0,T]) c \h\ c H. (4) 



and 



where, if if : [0, T]" 



E5 BH - K (uf<E\\u\\^ + E\\D BH - K u\\^ mnl (5) 



f d^R d^R 

IMIf«|®|*| = J \<f(u,v)\ \<p(u',v')\ -^^{u,u')-^-^{v,v')dudvdu' dv' . (6) 

We will use the following formulas of the Malliavin calculus: the integration by parts 

F5 BH - K (u) = 5 BH - K {Fu) + (D BH - K F,u) H (7) 
for any u G Dom(5 BHK ), F G D 1,2 such that G L 2 (Q;7i); and the chain rule 

n 

^•\(F) = £«F)I^V (8) 

i=\ 

if <p : M n — > R is continuously differentiable with bounded partial derivatives and i 7 = 
(.F 1 , . . . , .F m ) is a random vector with components in B 1 ' 2 . 

By the duality between F) bH K and $ bH ' k we obtain the following result for the con- 
vergence of divergence integrals: if u n G Dom(5 B ' ) for every n, u n — > u in L 2 (Q;?i) and 

i BH,K (ti n ) -> G in L 2 (ft) then 

« G Dom(5 B " ,K ) and 5 B "' K («) = G. (9) 
It is also possible to introduce multiple integrals I n (f n ), f G 7i® n with respect to B H,K . 

Let 

F = J2Wn) (10) 

n>0 

where for every n > 0, /„ G "H® n are symmetric functions. Let L be the Ornstein-Uhlenbeck 
operator 

LF = -Y^nlnifn) 

n>0 

if F is given by flTCTp . 

For p > 1 and a G M we introduce the Sobolev-Watanabe space O a ' p as the closure of 
the set of polynomial random variables with respect to the norm 



\F\L p = \\(I-L)2 



\LP{tt) 
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where / represents the identity. In this way, a random variable F as in (fTUj) belongs to B a > 2 if 
and only if 

l 2 (n) < °°- 

n>0 

Note that the Malliavin derivative operator acts on multiple integral as follows 

oo 

D? H - K F = Y j nI n -i{fn{;t)), t E [0,T]. 

n=l 

The operator Y) bH ' k is continuous from D a " 1,p into D a ' p (7i) . The adjoint of D bH ' k is denoted 
by S B ' and is called the divergence (or Skorohod) integral. It is a continuous operator from 
pa,p j n ^ j])0,p_ p or adapted integrands, the divergence integral coincides to the classical 
Ito integral. We will use the notation 

6 BH ' K (u)= ( T u s 5B s . 
Jo 

Recall that if u is a stochastic process having the chaotic decomposition 

Us = ^2 I n (/ n (-,s)) 
n>0 

where f n (',s) S 7^® n for every s, and it is symmetric in the first n variables, then its Skorohod 
integral is given by 

l-T 

I u s dB s = y^I n+1 (f n ) 

J ° n>0 

where f n denotes the symmetrization of f n with respect to all n + 1 variables. 

3 Tanaka formula for unidimensional bifractional Brownian mo- 
tion 

This paragraph is consecrated to the proof of Ito formula and Tanaka formula for the one- 
dimensional bifractional Brownian motion with 2HK > 1. Note that the Ito formula has been 
already proved in [15]; here we propose a different approach based on the Taylor expansion 
which be also used in the multidimensional settings. 

We start by the following technical lemma. 
Lemma 1 Let us consider the following function on [l,oo) 

Hy) = y 2HK + (y- lf HK - ^ (y 2H + {y- lf H f . 
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where H G (0, 1) and K <E (0, 1). Then, 

h(y) converges to as y goes to oo. (11) 
Moreover if 2HK = 1 we obtain that 

lim yh(y) = -{l-2H). (12) 

y— >+oc 4 

Proof: Let y = —, hence 

e 



e I e 2HK 



H(l-r-^(l + (l-efY 



Using Taylor's expansion, as e close to 0, we obtain 



h (}) = psj? ( h2r ( k ~ l > 2 + °( £2 )) ( 13 ) 



Thus 



lim h(y) = lim hil/e) = 0. 

For the case 2HK = 1 we replace in (113[> . we have 

1 fl\ 1, , 1 



-^-j= i( l-2 H) + ^) 
Thus (fl~2|) is satisfied. Which completes the proof. ■ 
Theorem 1 Let f be a function of class C 2 on M.. Suppose that 2HK > 1, £/ien 

/ ( S f >*) = / (0) + jf f {Bf K ) 5B»> K + HK jT /" (Sf s 2 ^- 1 ^. (14) 

Proof: We first prove the case 2HK > 1. Let us fix i > and let be ir := {tj = = 0, n} 
a partition of [0,t]. By the localization argument and the fact that the process B H > is 
continuous, We can assume that f has compact support, and so /, /' and /" are bounded. 
Using Taylor expansion, we have 

n i n o 

/ K' A ') = /(o) + E f K") - B t") + \Y,f «* - 

i=i j=i 
:= /(0) + /" + J n . (15) 
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where Bf K = B** + 0, - B^) with 6j is a r.v in (0, 1). 

Since B H,K is a quasi-helix (see [12]), we can bound the term J n as follows: 



E\J n \ 2 < nC/^ E { B t 3 ,K : -- B 5? V 



< 2- 2 ^nC^|t,-t J _ 1 | 4 ^ 

< 2 nCt svtp\tj—tj-i\ = 2 L 

3=1 

where C a constant depends de /". Then 

J" — ► in L 2 (0). 

n— >oo 

On the other hand, we apply ([7]) and we get 



4HK 



n AHK-2 n^oo 



(16) 



* B ' (^-i ) 1 (*i-i.*j](0 + X^" ( S *A ) ( 1 (o,t i -i]» 1 (t i -i,t J -]) 



J =1 



Next 



3=1 
n 

E 



2^ 



(( + 2H,+2H\K r + + \2HK\ + 2HK 
{{tj +t j _ 1 ) -(tj-tj-!) )-tj-l 



We denote by 



A, := HK / s 



ds = -t 



2HK 



To prove that Jg converges as n — > oo, it suffices to 

show that 

2\ 1/2 



C n :=\E 



0. 
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By Minkowski inequality, we have 



.7=1 V 



ii + 2H , ,2H \K /. , \2HK\ .2HK 



2<* 



n 



i=i 

= Ct™ K /2[Cl + Cl] 

where C is a generic constant. 

Since 2^ > 1 then C 2 n := ^k^nh^r 



1 - 2 1 

2lHf X/ I^C?)I + 2K n 2HK-l 



r 1 — 

■ n 2HK 



wk h( j) 



0. According to (jllf) . we obtain 



< 



3=1 



fflHK 1 rl __ >oc 



0. 



Thus 



We show now that 



/« _^ FA" C f" (B^ K )s 2HK - l dsmL 2 (n) 

rwoc Jo 



vi=i 



/' (B H «*)l (0>t] (.) in L*(H;H). 



Indeed, using the quasi-helix property of B H > , we obtain 



n 



3=1 

E tr r i/ , (<j)-/ , wii/ , (^f)-/ , (^i^(^)ducfo 

j 1=1 tj— 1 ■'*!—! 



— ^ /■*' ^i? 



< 2^(sup|r(x)|) 2 sup It,-^! 2 ^^ / / 

xeM i=l,...,n j i_-^Jtj—i Jti—x 

= 2 1 - K (sup\f"(x)\) 2 sup ^-^^(M) 
xeM i=l,...,n 



dudv 




(u, v)dudv 



This with JTH]) and (US} implies that ij» converges to 5 bH ' K (/' (5 H > A ') l (o t] (.)) in L 2 (ft). 
Therefore p^j) is established due to ©. ■ 



The proof of the case 2HK = 1 is based on a preliminary result concerning the 
quadratic variation of the bifractional Brownian motion. It was proved in [19] using the 
stochastic calculus via regularization. 
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Lemma 2 Suppose that 2HK = 1, then 



j'=i 



Proof: A straightforward calculation shows that, 



To obtain the conclusion it suffices to show that 



2fc-i ■ 



lim E(V t n ) 2 = (^t) 2 . 

n^oo Z 



In fact we have, 



n 



Denote by 



» n (i,j)=E((B? i > K -B^)(B? j > 
It follows by linear regression that 



JD 



»„(i,j)^i + \/Sn(i,j)-(0 n {i,j)) 2 N 2 



where iV"i and -/V2 two independent normal random variables, 



W H + 3 ) -2\j-i\- (f H + (j - 1) 2 T + 1| 



2 x n 

-((• - ir + J="T + \3 - i + 1| + ((• - I) 2 " + (J - 1) 



and 



Hence 



J::{B H, K _ B HXy E ^H,K_ B H^ 



Hn{i,j) = 2(6»„(z,j)) 2 + S n (i,j) 
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For 1 < i < j, we define a function fj : (1, oo) — > R, by 

) 2 



fiix) = {{x-lf H + f H ) K -{{x-lf H + {j-lf H ) K 

2H , -2H\K , /_2H , -i\2H\K 



We compute 



( x -iyH j y (x-i) 2H 

X 2H + j2 tf X A'-l , 2H + (j _ 1)2 ifs K-l 



+ 



:= g(x - 1) -g(x) > 
Hence /j is increasing and positive, since the function 



is decreasing on (l,oo). This implies that for every 1 < i < j 

l^.^l = ^/iW<2^/iO-)<^C7)l 

and \6 n (i,i)\ = ~\h(i) + 2| for any i > 1. 
Thus 

i,j=l i<j i=l 

Combining this with (fT2j) . we obtain that ^i\'=i ®n{hj) 2 converges to as n — > oo. 
On the other hand, by (fT2"j) 



E = J E f^) + 2^=r) (^') + pW) - 

ij=l m=i 



t ^ 2 



2 x-i 

Consequently, -E(V t ?l ) 2 converges to ( 2 k-i ) 2 as n — > oo, and the conclusion follows. ■ 

Proof: [Proof of the Theorem Q] in the case 2HK = 1.] In this case we shall prove that 

i? — f /' (*f «*) in L 2 (0), (17) 
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and 



To prove (|18p . it is enough to establish that 



Appliquing the Minkowsky inequality and (|12|) , we obtain 



„ n 

£ 



1 (tf 



2 A' 



< 



2n 



3=1 3=1 



By using Lemma [2l we conclude that 

The rest of the proof is same as in the case 2HK > 1. 



(18) 



(19) 



Let us regard now the Tanaka formula. As in the case of the standard fractional 
Brownian motion, it will involve the so-called weighted local time (x E R, t E [0,T]) of 
B H ' K defined as the density of the occupation measure 

A E B(R) — ► 2HK f l A {B^ K )s 2HK - l ds. 
Jo 

Theorem 2 Let (^Bf' K ,t E [0,T]j be a bifractional Brownian motion with 2HK > 1. T/ien 
/or eac/i t E [0, T], i£l i/ie following formula holds 

x\ + f sign{B s - x)5B^' K + Lf. (20) 

JO 



5 



t 
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Proof: Let p e (y) = ^ £ e ^ e be the Gaussian kernel and put 



and 



By the Theorem [T] we have 



F' £ {z) = 2 f p £ (y)dy ~ 1, 

J — oo 

F £ (z) = f F' £ {y)dy. 
Jo 



F £ (B?> K - x) = F £ (-x) + J q F> {Bf > K - x) SBf* 

+HK f p £ (B^ K - x) s 2HK ~ l ds. (21) 
J o 

Using (a slightly adaptation of) Proposition 9 in [19j , one can prove that 

L% = lim 2HK f p s (B?' K - x)s 2HK - l ds in L 2 {n). (22) 
Jo 

and L% admits the following chaotic representation into multiple stochastic integrals (here I n 
represents the multiple integral with respect to the bifBm) 

where H n is the nth Hermite polynomial defined as 

H n (x) = ^J- e x > 2 4- (e~ x /2 ) for every n > 1. 
w n! dx nX ' 

We have F £ {x) — > \x\ as e — > and since F £ {x) < then by Lebesgue's dominated convergence 

theorem we obtain that F £ (B^ ,K — x) converges to \B^' K — x\ in L 2 (J7) as e — > 0. 

On the other hand, since < F £ (x) < 1 and F' £ {x) — ► sign{x) as e goes to the 

Lebesgue's dominated convergence theorem in L 2 (0 x [0,T] 2 ;P® JuJ| ('"j v )dudv) implies that 

i 7 ^ (B H ' K — x) converges to sign [B H ' K — x) in L 2 {Q;TL) as e goes to because 

E\\F' e (B H > K -x)- sign (B H > K - x) ||f w , 

= £ / T T |F; - x) - Sl gn (B»> K - x) \ \F> (B*> k - x) - Sl gn [B»> K - x) \ 

Jo Jo 

92R ( SA A 

x [u, v)dudv. 
ouov 

Consequently, from the above convergences and Q 

f K {Bf > K - x) 5Bf K -> f Si gn ( B f K - x) 5B^ K in L 2 (tt). 
Jo E ^ u Jo 

Then the conclusion follows. ■ 
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4 Tanaka formula for multidimensional bifractional Brownian 
motion 

Given two vectors H = {H x ,...,H d ) £ [0, l] d and K = (K u ..., K d ) £ (0, l} d , we introduce 
the (i-dimensional bifractional Brownian motion 

B H,K = (B Hl,Kl ,...,B Hd,Kd ) 

as a centered Gaussian vector whose component are independent one-dimensional bifractional 
Brownian motions. 

We extend the ltd formula to the multidensional case. 

Theorem 3 Let B H,K = (B Hl,Kl , B Hd,Kd ^) be a d- dimensional bifractional Brownian mo- 
tion, and let f be a function of class C 2 (M. d , R) . We assume that 2HiKi > 1 for any % = 1, ...n, 
then 

f (ijf = / (0) + £ jf g (Bf «jf + g tf^ jf ' g (Bf a ^-icfa. (24) 

Proof: Let us fix i > and a partition {i,- = = 0, ...,n} of [0,t]. As in above we may 
assume that f has compact support, and so /, /' and /" are bounded. Using Taylor expansion, 
we have 



n a af 



3=1 i=l 
n d 



d 2 f 



B H ' K ) (B^' Ki - B^ Ki ) (B^ Kl - B H » K > 



dXidxi V 3 ) V *i 7 V *i 

j=l i,j=i 



:= /(0)+/™ + J n . 

where -B^'^ = B^l^ + #j (^B^' K — B^'^^J , and 0j is a random variable in (0, 1). 
We show that J n converges to in L 2 (Q) as n — > oo. 



ei /'f < „/4 E £ | £ ^ (*D (<"*' - <:f ) - <i'f' 



rv2 n d / x 2 

O « \ -> \ -> „ / „ rr. K, „ w„ jr.- \ z 



j=l i,«=l 



i,i=l j'=l i,i=l 
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According to 0, we get 

r 



n d 



j=i i=i 



E 

i=i 



j =1 



dxi 



3=1 



+ £ fog ( B h- K i) ^(o^-ii' ife-i,* 



E[ Tn,i , T-Ti.-i 



1=1 

As the similar way in the above theorem, we obtain that for every i = l,...,d 

"'f!((Bf,*) s 2H ^- l ds in L 2 (ft). 



J 2 ' — > HiKi 

n—>oo 



o dxj 



We show that for every i £ {1, d} 



We set 



By inequality (0), we have 



dxi 



(B?> K )l m (s). 



E[5 B " '(u n > 1 )) < E\\u n ' l \\^ +E\\Du 



.n.i II 2 



where is the Hilbert space associated to B Hi > Ki and i?j its covariance function. 
For every r,s <t 



3=1 * * 

we remark that D r iil ,% and it™' 1 converge to zero as n — > oo for any r,s < i. Since the 
first and second partial derivatives of / are bounded, then by using the Lebesgue dominated 
convergence theorem and the expression of the norm \'H t \ <S> we obtain that 



in L 2 (n). 
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The proof is thus complete. 



Set 



One can easily generalize the above theorem to the case when the function / depends 
on time. 

Theorem 4 Let f £ C 1 - 2 ([0,T] x R d ,M.) andB H ' K = [B Hl ' Kl , ...,B Hd ' Kd ) be a d- dimensional 
bifBm with 2HiK{ > 1 for any i = 1, n. Then 

f(t,B?> K ) = f(0,0) + J^(s,B^)ds + ^^^(s,B^)5B^ 

+i H ^I^U i ^ B ^ s2HiKi ~ lds - (25) 

We consider twice of the kernel of the d— dimensional Newtonian potential 

f r(d/2-i) i if^>o 

u(z) = \ " M d - 2 1 - 3 

\ \log\z\ if d = 2. 

n?=i^^7 v / 21ot / 

where x = (a:i, ...,x d ) G I rf and < 7 := ±(2 - d) + + (d - 2)(HK)* - Y%=l H i K i with 
(HK)* = max{H 1 K 1 ,...,H d K d }. 

We shall prove the following Tanaka formula. It will involve a multidimensional 
weighted local time which is an extension of the one-dimensional local time given by (|23p . 
Note for any dimension d > 2 the local time is not a random variable anymore and it is a 
distribution in the Watanabe's sense. 

Theorem 5 Let U as above and let B H,K = (B Hl ' Kl , B Hd,K A be a d- dimensional bifBm 
with 2HiKi > 1 for any i = 1, ...d. Then the following formula holds in the Watanabe space 

^2 _1 f or an V 01 < 2(HK)* ~ 

U(t, B?' K ) = £7(0, 0) + f d s U(s, B?> K )ds + E f dU{s f"' K) 5B^ + L% x) (27) 

J0 ~l Jo ° x i 

where the generalized weighted local time L s (t,x) is defined as 

L e (t,x)= y /*TT f^fr) H n J^L=)li(l® n \)s e ds. 

n=(ni,...,n d ) 1=1 ° 
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Proof: We regularize the function U by standard convolution. Put U £ = p d * U, with pf is 
the Gaussian kernel on R d given by 

d d l 2 

p d £ {x) = T[pe{xi) = TT -=e-fe, Vx = (xi, ...,a; d ) G R d . 
Using the above Ito formula we have 

Us(t,B?> K ) = U £ (0,0) + fi^(s,B?*)ds + £j^(s,B?*)8B?<" 



j7 e (0,0)+/f(t) + /i(t). 



On the other hand, if V(z) = ?7(ai£i, a^) and V E = * V we have 

1^ 1 d 2 V £ , , d , , 
2L^9jW =P e (oi«i,-,ad«d)- 



Hence 



m) = — d 1— — - fp d £ ( Cl (s)(B^ - Xl ), c d (s)(B^ K * - x d )) s 9 ds 

llj=l y ZHjKj Jo 

where Cj(s) = ^-!-j== for every i = 1, ...,c£. The next step is to find the chaotic expansion of 
the last term If. By Strook formula, we have 

oo 
n=0 n - 

and 

ED n p e ( Cl ( S )(B^ - Xl )) = «(s) n EpW {ci( 8 )(B?» K ' - x t )) lg; s] (.) 



2H i K i , e 



c-faVn'f & +s 2^x, r »/2 / v VV <*(«) / 

llj l Cj ( s ) 2+ j W ^+^ l ^ J Cj(s) n+1 IM^ 

n ' / £ , a 2HiKi\-nl2 n / \tt | ^ I 1®" M 
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Consequently 

pi { Cl {s){Bf^ - Xl ), ...M*)(B?**' ~ x d )) = e n 



CAS) 
n=(m,...,n d )eN d «=l 



and that 

n=(n lv ..,n d )eN d U i=l S 

rf d P a 2H i X i+ «_(a! i ) / 

V / TT ie^_ Hn x * r n .(i® n \)s e ds 



This term (in fact, slightly modified) appeared in some other papers such as Proposition 12 
in [10], or in |22j . Using standard arguments we obtain that the last term converges in ©2 to 
L e (t,x) as £ goes to 0, with a < 2 {hk)* ~ ^/^" 

The rest of the proof is to show that the following convergences are holds: 
For every i = 1, d 

f d t U £ (s, Bf C d t U (s, B*> K ) SB?***. (28) 
Jo Jo 

f 't ma ft 



[ d s U £ {s,B?' K )ds f d s U(s,B^' K )ds (29) 
Jo Jo 



and 



U&rf^-Z+UfrB"*"). (30) 

We start with the convergence (f28l) . Fix i G {l,...,d}, we note gl(s,z) = diU E (s,z). By the 
formal relation (<5 is the Dirac distribution) 



f(y)5{x - y)dy = f(x) 
we can write (this is true in the sense of Watanabe distributions) 

gi(s, B?> K ) = I gi(s, y)S {B?' K - y) dy. 
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Furthermore (see jlOj but it can be also derived from a general formula in |16| ^ 



i=l 
d 



i=l \ n >0 



n=(m,...,n d ) 

whereas) = Ri(s,s) = s 2H ^, A n (s,y) = Uti J^^PR^^HuAr^) and J n(lg| 



L [0,s] 



Ui=i p n , (^jpjj for ever y n = ( n i> -i"*)- 

Hence 



ri=(ni,...,n d ) 



gi(s,y)A n (s,y)dyI n (lf Qi 

®l' 
[o. 



n.= (ni,...,n d ) 

and using the chaotic form of the divergence integral 

Jf(i) := f gl{s,B»> K )6B^ K < 
Jo 



n=(m,...,n d ) 



n;+l 



^( S )if ;; ] ( S1 ,... )Sn ji M ( S )n^ (ig; 

2=1 



n=(ni,...,n d ) 



+ 1 



ni>0 



E /2( S l>-> S ni> S ) 
n=(ni,...,ni,...,n d ) 



where the superscript (s) denoted the symmetrization with respect to s±, . . . , s ni , s, and 

n.i+l _, d 
~ + 



Observe here that, since the components of the vector B H,K are independent, the term 
nj=i lnj (^-f) Si]) * s y i ewe d as a deterministic function for the integral l % n . . The convergence 
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is satisfied if the conditions i) and ii) of Lemma 3 in [TU] hold. It is easy to verify the 
condition i), we will prove only the condition ii). 
Fixing i € {1, ■■■,(£}, we can write, 



\\Jm\\l-l,2 < ^(m + ir- 1 £ (n i + l)\E 

m>l \n\=ni+...+nd=m— 1 

|n|=ni+...+nd=m— 1 
n;+l 



m>l 



f 



/ E 7 l n2 l^(g/)ll^(^)|l[o,t](^)l[o,i](^) 

[0,T]«i+ 1 7[0 i T]" l + 1 ^ («i + 1) 



3=1 



x -^-^ r (s q ,r q )dr 1 ....dr ni+1 dsi....ds rit+ i 

q=l 1 1 

Since a < then (m + 2) a_1 < (m + and (rtj + 1) < (m + 1). This implies that 



\\jf(t)\\l-i,2 < E( m+1 ) a E ("*)' 

|n|=ni+...+n d =m 



m>0 



(1 



[0,T] 2 ./[0,T] 2 



("i + 1) 

^• i ( Sl )||^ i (r 2 )|l M ( S i)l [0li ](r 2 )i? i ( S i,r 2 )^- 1 



(? 2 i? <9 2 i? 

]J njl-Rj^x, r 2 ) nj l[ , Sl ] (g2)l[o,r 2 ] ( r i) d j- Oi, r i) j^j^ (s 2 ,r 2 )ds 1 ds 2 dridr 2 
3& 1 1 



+ 



1 



cfe 2 dr 2 

|B^( S i)||^ < (n)|l [0 , t] ( a i)l[o, t ](ri)i2 j (fl 1 ,r 1 ) 



rij-i 



(Tlj + 1) 7[ ,T] 2 i[0,T] 2 

d d 2 R d 2 R 

^Qn j !^ j (si,ri) n n [0iS1 ](g 2 )l[o,ri]( 7 '2) ds ^ ^ (s 2 ,r 2 )ds 1 ds 2 dr 1 dr 2 

3=1 
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By integration we obtain 



\\jf(t)\\i-x,2 < ^2(m+ir J2 ( n *) ! 

\n\=n 1 +...+n d =m 

\B^( Sl )\\B^(r 2 )\R t ( Sl ,r 2 ) 



m>0 



(1 



1 



{rii + 1) ' 



n,;-l 



[OA 2 

d BR- dR- 

^n j \R j (si,r 2 ) nj -^{si,r 2 )-^{si,r 2 )dsidr2 

j¥* 1 



+ 



1 



(ra; + 1) J [0>t] i 
d 



{B^is^WB^ir^R^nr 
d 2 Ri 



TT njlRjfa, n) n ' — — — (si, ri)daidri 
XA dsidri 



3=1 



We have for any 1/4 < (3 < 1/2 



gl(s,y)A n (s,y)dy 



gi(s,y)l[H nj 



Vi 



1 e 



-/3y 2 j »| 

p £ d (*)<fe/ di U{s,y-z)Y[H nj [^=)- 7 =^=^ l d y . 



3=1 



Since (see Lemma 11 in [10J) 



sup n |iJ n . (^) \e~ p 'i <C\[ 



z£R c 



3=1 



J=1 ynj\{rij V 1) 12 
and for every (s,z) G (0, T] x R d 

^CsiC 1 -^ f( Zl - Xl )s' HlK1 ,...,(z d -x d )s- HdKd 



l-d 
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Then, for any s G (0, T] 

\B e /(s)\ < 



j=i \Jn,j\(nj V 1) 12 s 1 1 3 JR d 



/ s ' 2 

JR d 



\s HK (y — z — x) 



< Csr-ifl-—^ —^7rf PeWz 

fJi ^~\{rij VI)" s n ^ K ^ J Rd 



e y 2 



a-P)\s- HK yf 



\(y- z- x) 



-dy 



where s HK y := (s HlKl yi, s HdKd yd)- 

Let r\ a positive constant such that, for every s G (0, T], j G {1, . . . , d} we have s~ 2H i K i > rj. 



Combining this with for any a, b G R, a 2 > ^(a — 6) 2 — 6 2 , we obtain 



-(|-/3)| S ^| 2 



e v 2 



Rd e V |(y-«-x)| d - 



-dy < 



e /v 2 



^(|-/3)|y| 2 



e -l(i-/3)|2/-(^)| 2 

d=T d f 



JM d 

" L M ) L l(v -(* + *))! 

< ^(l-^NI 2 /" ^( z ) e 2^|-/3)k| 2 dz 



Hz 



R d y/2TT 



i_am^i2 / e 



-/?M 2 



R d \/2vr 



-civ < oo 



since 2-qe < 1 when e close to 0. 
Thus 



\B^(s)\<Cs^f[ 



where C is a constant depending only on (i, -ff, K, T, x and /3. 

Using this inequality a 2 + b 2 < 2ab, for every a, 6 G R+, we conclude that there exist a constant 
C(i?, K) positive, such that for every i = 1, . . . , d 



dRi dRi 

{r,s)— — (r, s) 



dr ds 
d 2 Ri 



drds 



(r,s) 



< C(H,K)(rs) 2H ^-\ 
< CiH^K^rs) 11 ^- 1 
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and 

Ri{r,s) 

It follows by anterior inequalities that 



< C(H,K). 



ll^(t)llLi, 2 < C5> + m) a E tl~, ^TIO 

m>0 |n|=ni + ...+n d =mi=l ( n j V 1) 6 



[o,t]2 (rs)^- 1 )^ 



iKi 



n (,, s)i 



drds. 



We use the selfsimilarity of the covariance kernel R(r,s) = R(l, f.)r 2HK and the change of 
variables r/s = z in the integral respect to dz to obtain 



H«/fWllLi, 2 < Cr£(l + m) a E II ^ET 

■m>0 |n|=m+...+n ti =m j=l ( n i V !) 6 

-1 / n /i . \ \ ru-l d 



X / r 2(7-l+^^) dr / 



n 



Since for each i G {1, . . . , d}, 7 — | + i/j-ftTj > —1, 2(7 — 1 + H-iKi) > —1 and from Lemma 12 
and the proof of the Proposition 12 in [TO], we obtain that 

\\Jf(t)\\l-l,2 < C^ + mTm-^T E II- ~7mET 

m>0 |n|=ni + ...+n d =mi=l ( n j V !) 6 



m>0 



and since is finite if and only if a < 2 {hk)* ~ I ' 

On the other hand, for every (s, z) S (0, T] x M d we have 



d.U(s,z)\ <Cs^ +t) Uz 1 -x 1 )s- H ^ Kl ,...,{z d -x d )s- n * 



-H d Kd 



and 



U(s,z)\<Cs^ 2 -^+ 9 Uz x -x x )s- H ' K \...,{z d -x d )s 



H d Kd 



2-d 



2-d 



and this inequalities imply as in [21] the convergences (f29j) and ([30]) in mathbbD^ ■ 
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